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Abstract 

An approximate strategy for studying the evolution of binary systems of extended objects is 
introduced. The stars are assumed to be polytropic ellipsoids. The surfaces of constant density 
maintain their ellipsoidal shape during the time evolution. The equations of hydrodynamics then 
reduce to a system of ordinary differential equations for the internal velocities, the principal axes 
of the stars and the orbital parameters. The equations of motion are given within Lagrangian and 
Hamiltonian formalism. The special case when both stars are axially symmetric fluid configurations 
is considered. Leading order gravitational radiation reaction is incorporated, where the quasi-static 
approximation is applied to the internal degrees of freedom of the stars. The influence of the stellar 
parameters, in particular the influence of the polytropic index n, on the leading order gravitational 
waveforms is studied. 
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1 Introduction 

The search for gravitational waves is one of the most challenging projects of 21st century physics. In- 
spiraling compact binaries are among the most promising sources for gravitational waves that could be 
detected by gravitational waves observatories, such as LIGO, VIRGO and GEO600. During most of 
the inspiral time the gravitational waves emitted by a compact binary system are much to weak, the 
frequencies are much to small to be detectable by existing gravitational wave observatories. However, at 
orbital distances corresponding to the last 10 or 15 minutes of the inspiral time, the orbital frequency 
increases from around 1 Hz up to 1000 Hz, while in the same time the amplitude of the gravitational 
waves increases, thus improving the chance of actually detecting the weak signals. 

In order to extract a possible signal from the detector noise the gravitational waveforms emitted by 
inspiraling binary systems must be known in great detail, in particular during the last stages of the 
inspiral process before the final plunge. The frequency gap covered by contemporary gravitational wave 
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observatories ranges from 1 to 1000 Hz. In this regime post-Newtonian (pN) effects become important 
and have to be included into the analysis. For non spinning point-particle binaries an analytic solution is 
available up to the 3rd pN approximation [IJ. Incorporating the spin leads to enormeous complifications. 
Until now there exists a solution only for point-mass binaries with either two equal masses and arbitrary 
spins, or two arbitrary masses and only one spinning object 

However, while the components of black hole-black hole binaries can be treated as pointlike objects even- 
tually up to the innermost stable circular orbit (ISCO), yet another effect has to be incorporated into the 
analysis for all other binary systems: the internal structure of the components. Lai and Wiseman argue 
that during most of the final inspiral process even the neutron stars can be treated as pointlike objects 
13] • This is true if one considers only a few orbits. If, however, one is interested in the long-term evolution 
of the system, the tidal interaction of the neutron stars will lead to a phase shift in the gravitational 
waveforms, which is not neglegible anymore. 

In this paper we shall focus on the problem of including tidal interaction into the equations of motion 
and we discuss the influence of the internal stellar dynamics on the binary's dynamics as well as on the 
gravitational waveforms in leading order approximation. Post-Newtonian effects will be included in a 
forthcoming paper. 

A detailed description of the tidal interaction and the corresponding changes in the gravitational wave- 
forms requires the application of three-dimensional numerical hydrodynamics. This is well beyond the 
scope of this paper. Instead we shall consider an approximative formulation of the problem. This strat- 
egy allows an analytic formulation of the equations of motion. Moreover, we are able to run long-term 
calculations of the system's dynamics with sufficient accuracy. There are basically two approaches for 
an approximate description of the tidal interaction. The first one is based on a linear adiabatic theory, 
describing small derivations from equilibrium. Assuming an adiabatic equation of state, the linearized 
equations of hydrodynamics have been derived by Ledoux and Walraven 0] and Dziembowski 0- Press 
and Teukolsky j^j developed a technique for calculating normal modes of non-radial stellar oscillations 
in the presence of tidal forces. Using the formulation of Dziembowski, Kokkotas and Schafer studied the 
leading order dynamics and gravitational waveforms of an NS-NS binary system, in particular taking 
into account the influence of tidal interaction on the actual waveforms 0. Resonant tidal excitations of 
a binary neutron star with polytropic equations of state have been investigated by Ho in leading order 
approximation [Sj. 

However, there are scenarios when the linear adiabatic theory is not applicable. In particular the theory 
fails when the stellar oscillations cannot be considered as small derivations from equilibrium. Thus far 
only little investigation has been undertaken in this direction. For a rotating, oscillating disk of dust 
there exists an analytic solution up to first post-Newtonian order jS] . A system of an rotating, osillating 
dusty disk and a point-mass object has been studied by the author in a recent paper ^Hj- A so-called 
'affine stellar model' for polytropic, ellipsoidal configurations has been developed by Carter and Luminet 
|llj . Here the density contours in the star are assumed to form homologeous ellipsoids. The model 
allows for the discussion of arbitrary large amplitudes of the tidally generated oscillations, but fails to 
give accurate results for small oscillations. The excited mode roughly corresponds to the /-mode of the 
oscillations (see also Later on Lai and Shapiro extended this model to allow for the investigation 

of Riemann-S binaries |13|. In this paper we shall follow the approach of Lai and Shapiro, pointing out 
some discrepancies in the authors' formulation of the equations of motion as well as for the gravitational 
waveforms. 

The paper is organized as follows: In section 2 we briefly recapitulate the Lagrangian formulation of the 
dynamics of Riemann-S binaries derived by Lai and Shapiro in a series of papers ^3]> |14| . |15j . In this 
model the stars are considered as rotating and oscillating triaxial ellipsoids with a polytropic equation 
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of state. To keep the analysis as simple as possible, the tidal interaction potential is truncated, only the 
leading order (i.e. quadrupole) tidal interaction is taken into account. For later reasons we shall also 
present the Hamiltonian formulation of this problem. In particular we shall specialize to axially symmet- 
ric polytropic stars 1 . Section 3 is devoted to the derivation of the leading order gravitational reaction 
terms appearing in the dynamical equations. We shall point out some discrepancies between our approach 
and the approach taken by Lai and Shapiro. The leading order gravitational waveforms are derived in 
section 4. In particular we shall apply the quasi-static approximation to the internal stellar dynamics. 
The corresponding equations are derived for general triaxial Riemann-S binaries, but for simplicity we 
restrict ourselves to axially symmetric configurations in the numerical calculations. The numerical results 
and some possible physical applications of our model are discussed in section 5. In particular we study 
the influence of the polytropic index and hence of the equation of state on the dynamics as well as on 
the leading order gravitational waveforms. We shall show that the effect of the internal stellar structure 
on the long term evolution of the binary is quite significant. We argue that the exact knowledge of the 
polytropic equation of state is essential in order to calculate the leading order gravitational waveforms 
emitted by the binary system correctly. 

2 The Newtonian dynamics of Riemann-S binaries 
2.1 Lagrangian formalism 

On Newtonian order the Lagrangian equations of motion for compressible Riemann-S binaries have been 
derived by Lai and Shapiro in several papers (see |13j . |14| . and we shall therefore only shortly 

recapitulate their strategy. We then turn our attention to the Hamiltonian formulation of the problem 
which will be more suitable for numerical calculations. The case of axially symmetric, polytropic stars 
is of particular interest for us, but requires some special care in order not to overcount the degrees of 
freedom. 

Consider a binary consisting of two ellipsoidal fluid configurations M and M', asssuming the stellar 
matter to obey a polytropic equation of state, 

P = Kp 1+1 ' , \ P' =K'p a+1 l n '. (1) 

Here n and n' are the polytropic indices, while K and K' represent constants which are determined by 
the equilibrium radii of non-rotating, spherically symmetric polytropes of the same mass. A Riemann-S 
ellipsoid is characterized by the angular velocity Q — Qe z of the ellipsoidal figure around the principal 
axes and the internal vorticity £ =£e 2 . The model is based on two assumptions: First, we shall require 
the vorticity to be uniform, and second, the surfaces of constant density inside the stars are assumed 
to form self-similar ellipsoids. In other words, each star is described by only five degrees of freedom: 
the three semi- major axes 01,02,03 and a^a^a^ respectively, and the two angles tp and 7, which are 
introduced in figure 1. 

Consider now an isolated star M. If the assumptions given above are to be fullfilled, it is immediately 
clear that the fluid velocity in the body fixed system, which rotates with angular velocity = Qe^ 
(fi = 7) with respect to an inertial system centered at M, obeys the following ansatz: 

u c = Qix 2 ei + Q2Xie 2 . (2) 

1 Just setting ai = 02 in the equations of motions gives rise to divergent terms. 
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Here Q±, Qi are constants and x±, x 2 denote coordinates in the body fixed frame. The second requirement 
yields 



Ql — C — A, Q 2 — 2 2 2 C = A) 

a\ + 0,2 a x + a\ 



where C is the vorticity in the body fixed frame and A = ip is the angular velocity of the internal fluidal 
motion. On the other hand, observing the stellar matter in an inertial system centered at M it's velocity 
is given by 

ujs = u c + fi A r. (3) 

The rotational energy of the star reads 

Trot = ^ J pu IS ■ U IS dV 

k M 9 

= ^-{a\ + al){K 2 + n 2 )- Z -K n Ma 1 a 2 Kn, (4) 

where K n < 1 is a constant, which can be obtained from the Lane-Emden equation (see Appendix). It 
can be easily seen that the inertial tensor iy = J pXiXjdV in the body fixed system takes a very simple 
form, 



5 QllJ ' K "-"3 ^l^il 



Iij = / pXiXjdV = ——aidij, K n := - u 4 — . (5) 



In particular, for a\ = 02 the rotational energy reduces to 

K n M 



T rot = ^—{n-KYa{. (6) 
5 

The total kinetic energy of star M is given by 

T s = ^f(aj + al)(n 2 + A 2 ) - ^Ma^AQ + ^(d? + d| + a 2 ), (7) 
and the Lagrangian L s of the star reads 

L S =T S -U- W. (8) 

The gravitational self-energy W , and [7, the internal energy of the star, have been computed by 13 . 
They are given by 

U = J n^dm = kxKMpl /n , (9) 

n(n+l) 



where k\ = ™_ n is constant, and 
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Here R — (aia 2 a 3 ) 1 / 3 is the mean radius of the ellipsoid and J is given by J = a\A\ + a 2 A 2 + afA 3 . 
The coefficients Ai, defined by Chandrasekhar ^Sj, are given in the Appendix. 

Now let us come back to the binary system. The tidal interaction potential is clearly dominated by 
the quadrupole interaction, which will be the only terms to be included into our calculations. Thus the 
orbital Lagrangian reads 

L rb = ^r 2 + ^4> 2 -W int , (11) 

where p — MM' /(M + M') is the reduced mass, and the interaction potential Wi n t is 

GMM' GMM'n n r „ , 9/ , „ 

Wint = t-— ^ — — [a? 3 cos 2 a - 1) + a\ 3 sin 2 a - 1) - a 2 ] 

_ GM ^ K ' n [a? (3 cos 2 a' - 1) + a 2 2 (3sin 2 a' - 1) - a?] . (12) 

The meaning of the angle a = (j> — 7 can be read off from Fig. 1. It is now easy to write down the 
Lagrangian of a general Riemann-S binary according to 

p , 2 ^r 2 - 2 GMM' GMM'n n r 2 ^ o 2 2 . 2 2 , 

L=-r +—0 + — - — + — — [af(3cos a - 1) + a 2 (3sm a - 1) - a 3 )J 

+ [ a f( 3 cos 2 a' - 1) + a 2 2 (3 sin 2 a' - 1) - a 3 2 ] 

+ ^L( a \ + a 2 )(A 2 + f! 2 ) - -«„Af ai a 2 Ar! + ^(a 2 + a 2 + d 2 ) - k x KMp l J n 
10 5 10 

+ -^^f^ + ^K 2 + «2 2 )(«' 2 + a' 2 ) - 2 - K ' n M'aW^n' 

+ <§- (df + « 2 2 + d 3 2 ) - fc^'M'^"' + 5 ° n , (13) 

It is straightforward to derive the equations of motion governing the dynamics of the binary system. Let 
us first focus on the most general case when all three semi-major axes ai are different. With the central 
density p c being proportional to 1/ (0,10,2(13) and using 

dCLi di 

(see Appendix) the equations of motion derived from the Lagrangian l|13|) read 

GM' , , , , , , fhkxP c \ 1 3GM aiAi , , 

Sl - -W a ^ 3cos a - ^ + fll ^ + - 2a2Afi + \~£) a x ~ MT^t) W' (14) 

d 2 = ^ a2 (3sin 2 a - 1) + a 2 (A 2 + ^ 2 ) - 2 Ol A0 + f ^ 1 - (15) 



wt„ p c / a 2 k„(1 - §) 2R 3 



^_ GM' ^ | / 5fci P e \ 1 3GM a 3 A 3 
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i> = A = 



7 = 



Ij - u, 2 

1 



3GM' 



r 

3GM 



<,j — u, 2 

3 GM r 



a\a,2 sin 2a — 2(aidi — a 2 a 2 )A — 2(dia2 — aid2)fJ 

/i i 

(a 2 + a 2 ) sin 2a — 2(aidi — a 2 a 2 )£l — 2(axa 2 — dia2)A 

2 



, -1 K,Ja?(3cos 2 a - 1) + al(3sin 2 a - 1) - all 
r 1U r L 

+ «' n [a? (3 cos 2 a' - 1 ) + a 2 2 (3 sin 2 a' - 1 ) - a%] } , 

-< (a 1 — a 2 )K n sin 2a + [a{ — a 2 )K n sm2a >, 



10 r 5 



(17) 
(18) 

(19) 
(20) 



with M = M + M' . The corresponding equations for a' i} ip 1 and 7' are obtained by replacing unprimed 
variables by primed ones. Note that these equations apply whenever a\ ^ a 2 . If the two semi-major 
axes a\ and a 2 are equal, a x = a 2 , some special care is needed. In particular, Eqs. i|17fl an d l|18|l are ill 
defined, indicating that ip an d 7 are n °t independent variables. In fact, setting ax = a 2 , the Lagrangian 
(|13|l simplifies to 



H , 2 ^r 2 - 2 GMM' GMM' 
L= 2 r + — * + ^— + ^0rT 



I 1 2 \ 1 / / /2 /2\1 1 K n,M 2 • 2 
.«n(Ol - <4) + K ni a X - a 3 )J + —5— Ol/ 3 



1 K 'nM' II 'all 1 K nM 2 . -2\ , K 'n^' fr> ■ 12 . ■ I2\ i j\ r V™ 

+ — a i P + ^7j~( 2a i + a s) + ~~ [7r~( 2a i + a 3 ) - kiKMpJ 



5-n 2i? 3 



5-ra' 2i?' 3 



(21) 



where we introduced a new variable /3 := 7 — 1/). The number of degrees of freedom of the binary system is 
thus reduced from 12 to 8, resulting in an enormeous simplification of the numerical calculations. Using 
Eq. (|94() , the equations of motion for binary systems with axially symmetric stars become 



ax 



«3 



GM' 



2r 3 
GM' 



ax + ax$ + 



5fci P c \ 1 



3GM axAx 



r = 


rcf> 2 




r 






4> = 






r 


P = 


-2^/3. 




ax 



a 3 



GM 



K n n p c ) ax K n (l - |) 2i? 3 
5kx P c \ 1 3GM a 3 A 3 



nK n p c J a 3 K n (l - |) 2R 3 ' 
3 GM r 



c„(a 2 -a 2 ) + <(ai 2 -4 2 )}, 



(22) 
(23) 
(24) 
(25) 
(26) 



Later on we shall study the gravitational waves emitted by such particular systems in detail, but first we 
shall derive the Hamiltonian formalism for general Riemann-S binaries. 



2.2 Riemann-S binary systems in Hamiltonian formalism 

In the previous section we have revisited the derivation of the Lagrangian equations of motion for a 
Riemann-S binary. For the following analysis it is, however, more suitable to apply the Hamiltonian 



(i 




Figure 1: The figure shows the coordinate systems used throughout this paper. The a;j-axes of the 
comoving system £ is chosen along the axes between the centers of M and M' and includes with the 
X— axes of the inertial system the angle cf>. Dotted lines represent the axes of the body-fixed system. 
COM denotes the center of mass. 



formulation of the problem. The generalized momenta pi — dL / dqi derived from the Lagrangian are 
given by 



Pr 

Pa, 



[if, p^ = fir 2 (j), 
K n M . 



(27) 



K„M r 



Pi> 



(a\ + a 2 )tp — 2aid27 



K n M r 



Pj 



(a\ + a\)^ — 2a 1 a 2 ip 



For a\ ^ a 2 it is possible to invert these equations in order to express the generalized velocities in terms 
of generalized momenta: 



7 



Pr 




V± 


A* 


<}> = 




5 


r 2 


Ya\ 


K n M 


[(«?- 


4) 2 



p~f 



a\a 2 



Pa,, 



Pu- 



(28) 



[a\-alY 



K n M 



P^> 



a\a 2 



{a\-alY 



P- 



The corresponding equations for tp' and 7' are obtained by replacing unprimed variables by primed 
ones. For a\ = a 2 the generalized momenta p^ and p 7 are not independent variables, as can be easily 
seen from Eqs. (|27|> . This is reflected by the observation that the inversion problem Pi(qi) — > qi(pi) is 
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ill defined in this case. We shall turn to this later on, but for now we focus on the most general case, 
assuming a\ ^ a 2 . A Legendre transformation of the Lagrangian (|13l) leads to the Hamiltonian 



K 3 K 3 

Tp 2 + — Vp' 2 



a? + a 2 



2/z ' 2pr 2 ' 2n n M j-^^ ' 2K' n M' ^-^ ai ' 2n n M (aj - a%) 2 ^~< 



(p 2 7 + pI) 
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a l a 2 -J "-1 T- u-2 / /2 1 /2\ 1 ±u "l"2 / / 

K n M (a? - a 2 ) 2 P ^ + 2 <M' « 2 - a 2 2 ) 2 Vi ' ] K> n M> (a 12 - a' 2 ) 2 ^ 

3 GM 2 



a' 2 + a' 2 



k 1 KMp i 

GMpn n 
10r 3 

GMpn' n 
10r 3 



l/n 



2i? 3 



J + k' x K'M'p'}l n ' 



3 GM' 2 _ GA4p 



5-ti' 2i? 



/3 



a 2 (3 cos 2 a - 1) + a 2 (3 sin 2 a - 1) - a 2 



af (3 cos 2 a' - 1) + a 2 2 (3 sin 2 a' - 1) - a% 



(29) 



Once the Hamiltonian is given, the corresponding equations of motion 

dH dH 

opt oqi 

are easy to calculate using the relation for djjdai derived in the appendix. With the equations for the 
qi already given in Eqs. (|28|l the remaining equations read 



Pai 



Pa 2 



a\(ia\ + a?) . 2 2 \ a 2 . + 3a? 

(a2 _ ai)3 <Pt+P*) + 2a2 W 



kiMP c 3GM 2 ai Ax 



GMp 

5r 3 



K n ai(3 cos 2 a — 1) 



nax p c 



5 - n 2R 3 ' 



a 2 (3al + a|) 
(af-al) 3 ^' 



(of - a 2 .) 



?tyP 7 



GA4p 

5r 3 



K ra a2(3 sin 2 a — 1) 



kxMP c 3GM 2 a 2 A 2 
na 2 p c 5 — n 2R 3 ' 
fciM P c 3GM 2 a 3 A 3 GAf/i 



na 3 p c 
2ty = 0, 

3 GMp 



5-7i 2i? 3 



5r 3 



p 7 



10 r 3 
GMp 



r" pr 
3 GMp< 



ftn(<Zi — a 2 ) sin 2a, 
3 



3 GMpK n 
10 r 3 



a 2 (3cos 2 a - 1) + a?.(3sm 2 a - 1) - at 



10 r 4 
3 GAfp 
10 r 3 



af (3 cos 2 a' - 1) + a 2 2 (3 sin 2 a' - 1) - a[ 



[n n (a\ - a 2 ) sin 2a + «4(af - a 2 2 ) sin 2a'] 



(30) 



(31) 

(32) 
(33) 
(34) 



(35) 
(36) 



Together with the corresponding equations for the primed variables this system of differential equations 
describes the dynamics of the binary completely. 
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Let us now assume a% = a 2l a[ = a' 2 . As we have already pointed out in this case the dynamical 
equations l|28|) and (|30|l - l|36() do not apply. The generalized momenta derived from Lagrangian (|21|l read 



Pai = -KnMax, 
5 

Pr = pr, 

which can be easily inverted, giving 

5 



Pa 3 = — ^«3, 

5 

pr 2 <j>, 



P0 = -KnMa^, 



(37) 



ai 



■P ai : 



03 



Pa 3 > 



(3 



5 P/3 

2n n M a\ ' 



2n n M 

. Pr ; P<p 

r = — , 4>= 

p pr* 

As before a Legendre transformation of the Lagrangian (|21|l leads to the corresponding Hamiltonian 



(38) 



H 



2n n M 



2 +Pa3 + 2af 



2k' M' 



n 



T>' 2 



2a? 



Pi , pI 



2p 2pr 2 



kiKMp 



<SGM^3_ /B , _ 3GM- _GMp_GMp 

5 - n 2B? +klK M Pc 5-n' 2i?' 3 r 10r 3 LMfll ^ + M&1 ° 3 jJ ' 

(39) 

Note that /3, /?' and </> are cyclic variables, i.e. the corresponding generalized momenta are constants. The 
Hamiltonian equations are then given by Eqs. (|38|1 and 



P ai 

Pa 3 
Pr 



5 P% GMpK n n k x MP c 3GM 2 a 1 A 1 
01 + 2 ^3-, 

nai p c 5 — 73 K J 



2K n M a\ 5r 3 

GMuK n k x MP c 3GM 2 a 3 A 3 
■ a 3 1 



5r 3 
GMp 



na 3 p c 
p\ _ 3 OM/i 
/xr 3 10 r 4 



5-77 2i? 3 



[«„(a 2 - a\) + 4(af - a* 2 )] , 



P<t> = Pp = 0, 

and, of course, the corresponding equations for p' a . and j3' . 



(40) 

(41) 

(42) 
(43) 



3 Leading order radiation reaction in Riemann-S binaries 

On the Newtonian level, tidally coupled binaries form a conservative system. However, according to 
the theory of General Relativity binary systems loose energy due to the emission of gravitational waves, 
the first nonvanishing dissipative terms appearing at 2.5 post- Newtonian approximation. The following 
section is devoted to the calculation of the leading order radiation reaction terms for Riemann-S binaries. 
This can be done by virtue of the Burke-Thorne radiation reaction potential $ rea c, given e.g. in |17| . 



" ac _ 5c 5 ah XaXb ' 



(44) 
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As before the x a denote coordinates in the corotating coordinate frame. Thus the computation of the 
radiation reaction potential reduces in principle to the calculation of the 5th time derivative of the 
STF mass quadrupole tensor in the corotating system. The mass quadrupole tensor being additive at 
Newtonian order, it is possible to consider orbital and stellar contributions separately. 



3.1 Time derivatives of the stellar mass quadrupole tensor 

According to Eq. I|44|l the leading order gravitational wave emission is governed by the time variations 
of the STF mass quadrupole tensor. Thus even an isolated, oscillating star represents a source of gravi- 
tational waves. In a binary system tidal coupling between stellar and orbital degrees of freedom leads to 
a change in the gravitational wave pattern. Though the orbital contribution in general clearly dominates 
the gravitational waveforms emitted by the binary, the stellar contributions cannot be neglected. 
Let us consider a coordinate transformation such that the origin of the corotating system is centered at 
M. The relation between the coordinates X a of an inertial system centered at M, the coordinates x a of 
the body- framed system and the coordinates x a can be read off from Fig. 1, 



where 



X a = T ab (~/)X b , X a = Tab(4>)Xb, 



cos (j) sin (j) 
T a b{(t>) = — sin0 cos0 
1 



(45) 



(46) 



As has already been mentioned, the STF mass quadrupole tensor takes a particular simple form in the 
body- fixed system, 



lab — lab — g<$afc7 c , 



lab = / pX a X b dV = 



a 2 a 5ab- 



(47) 



To calculate the time derivatives of the mass quadrupole tensor in the corotating system we could in a 
first step consider the time derivatives of the star's mass quadrupole tensor in the inertial system, the 
later one being related to l a b by the transformation 

and only in a second step apply the transformation into the corotating system according to 



i^ = T- aa {mM l % IS) - 



(48) 



This calculation is straightforward, but successively inserting the Newtonian equations of motion leads 
to rather complicated expressions. In our approach another strategy is more suitable. Following Lai and 
Shapiro (see also ^Hj) we combine the steps mentioned above according to 



d 5 



I<J =T aQ (0)T 5 ^)— 71(7)7^.(7)^ 



E 

5 m 

EE 

m— p— 



dt 5 
d 5 ~ 



dt 5 - 



E 

p=0 



m 
P 



d m ~ p + 
7 1 aa (<P)-^Z^ T li(l) 



dP 



5 
m 



m 
P 



r 



dt 5 ~ 



(49) 
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where 



The calculation of the matrices R?, is straightforward and given in the Appendix. However, the resulting 

(5) 

expressions for being rather complicated, we should imply further assumptions on the internal stellar 
motion. In particular, it is reasonable to consider all internal velocities and accelerations to be small, 
thus applying the quasi-static approximation to the stellar degrees of freedom. This strategy was already 
followed by Lai and Shapiro J^jj but they applied the quasi-static assumption to the orbital motion, too. 
While this is justified for circular orbits or if one is interested in a few cycles only, for elliptical orbits this 
leads to a growing phase error in the gravitational waveforms. Since we are interested in the long-term 
evolution of the system we shall apply the quasi-static approximation to the stars only. Neglecting all 
terms of order 0(di) and 0(£l), keeping only terms linear in hi and fl and using that \a,i\ <C \^ai\ , Eq. 
(I49|l simplifies to 



p=0 



/c< p 



i 



4 
p=0 



(50) 



Note that I a {, w in this approximation. After some algebra we end up with 
ig 5 = (In to) 



40ft 4 (A i 




cos 2a 
— sin 2a 




80fi d fi 



- sin 2a 
cos 2a 




- cos 2a 
sin 2a 




sin 2a 
cos 2a 




(51) 



It is remarkable that the components of I_g are nonvanishing only for a\ ^ a.2- In other words, the quasi- 
static approximation does not allow for gravitational wave emission of an isolated, axially symmetric, 
poly tropic star 2 . 



3.2 The orbital contribution to i 



(5) 
ab 



As has been outlined before, the quasi-static approximation does not apply to the orbital motion in 
general. To calculate lj^'° rb we compute l^'° rb in the inertial system by successively inserting the 
Newtonian equations of motion 



GM 



r<t> 2 



rd> 
-2 — . 



(52) 



More precisely, one should insert the equations of motion of the tidally coupled system (cf. Eqs. (|19|l 
and (|20H ). This would complicate the problem enormeously. However, these corrections, being of higher 

2 Without the quasi-approximation there exist nonvanishing contributions to I^j for axially symmetric polytropes, too. 
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order in 1/r, can easily be neglected. 

In the inertial system the nonvanishing components of l\° read 



I 



(orb) _ \±r 
11 _ 6 



(1 + 3 cos 20), 



I 



(orb) _ 
22 - c 



(1-3 cos 20), 



j(orb) 
*33 



fir 
IT 



I 



(orb) _ fir 
12 — o 



sin 2(f>. 
(53) 



Using Eqs. (|52(l to calculate in the inertial system and finally transfering to the corotating system 
centered at the center of mass according to 



we end up with 



j(5,orb) 
db 



T- aa (4>)T- b0 (w:> 



(5),(orb,IS) 



(54) 



I 



(B),orb 



8GM[i 



I 



(5) ,orb 



2GM\x r 



GM n o -,0 

4 —— + 3r 2 + 18r 2 2 
r 

GM 



T (5),orb _ 2GMji r 

33 - 3 r 4 

t(5), orb . GMfl 

hi -" 4 ^3— ' 



r 

GM 

r 

2GM 

T 



+ 6f 2 + 81r 2 6 2 



+ 6f 2 - 9r 2 2 



9r 2 - 6r 2 <i 2 



(55) 



3.3 Leading order gravitational radiation — general case 

In a recent paper we discussed the dynamics and gravitational wave emission of a binary system consist- 
ing of a rotating, oscillating dusty disk and a point mass |l()j . There we incorporated the leading order 
radiation reaction into the Hamiltonian equations by adding a dissipative part to the Hamiltonian. The 
Hamiltonian equations can then be applied in the usual way. Here we choose an alternative approach, 
leaving the Lagrangian or Hamiltonian unchanged but modifying the Euler-Lagrangian equations accord- 
ing to 



d dL _dL 
dt dqi dq 

The generalized dissipative forces T qi are calculated from the energy dissipation rate 



W = - / V- V$reacpdV (57) 

as T qi — . Let us consider the contribution Wm of star M to the energy dissipation rate of the binary. 
The velocity v of a fluid element of M can be separated according to v = u + u or f,, where 

/ a i « ( a 2 . „\ 0,1 0, 2 &3 

u = — A - \l x 2 ei + A -I- il xie 2 H xiei H x 2 e 2 H a; 3 e 3 



a,2 J \ a\ J a\ a 2 03 

is the velocity of a fluid element relative to the center of M and 

u orb = -r M ei - r M i>e 2 (58) 
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is the orbital velocity of the star's center of mass. As before the ej denote the unit vectors in the body- 
fixed coordinate system, while the e a represent unit vectors in the corotating system. As one can easily 
read off from Fig. 1, the coordinates are related according to 



X\ = X\ cos a + X2 sin a — rj./, 



X2 



-x\ sin a + X2 cos a, 



%3 = x 3 . 



(59) 



Inserting this into Eq. I|57|) and using J pxiXjdV — K "^ 1 afSij the contribution of the star M to the 
energy dissipation rate yields 



2G K n M r 



5c 5 



a\CL\ (iiJ cos 2 a + li 2 sin 2 a — iiJ sin 2a) + a2a 2 sin 2 a + tkJ cos 2 a + liJ sin 2a) 

1 



+ n(al - a 2 2 )(i^ cos2a + ~(lg> - sin2a) + a 3 ds 



2GM r 



5c 5 



• iW 1 2 ii d 



(60) 



The corresponding contribution of M 1 is obtained by replacing unprimed quantities by primed ones. Note 
that according to Fig. 1 the orbital velocity of M' is given by u' orb = r' M ei + r' M (j)B2- Adding up both 
contributions and using Mtm^m + M'r' M f' M = firr, the energy dissipation rate of the binary system 
reads 



w = 



2G K n M 1 



aidi (l W cos 2 a + I~ sin 2 a — I ^ sin 2a) + 0202 sin 2 a + 1^ cos 2 a -I- sin 2a) 



r( 5 ). 



( 5 ). 



r( 5 ). 



l(5), 



+ 0(a? - a 2 ) (l^ J cos 2a + -(1^ - sin 2a) + I^ ; a 3 d 3 



2G K' n M' 
5? 5 



c^di (1^ cos 2 a' + I^ sin 2 a' - f § sin 2a') + a 2 d 2 (l^ sin 2 a' + cos 2 a' + f § sin 2a') 



n(af - a 2 2 )(l© cos2a' + 1(1® - lg } ) sin2a') + [g a' 3 a' 3 



5 C 5 v 11 ~ ^ 12 > 
It is now straightforward to derive the generalized forces 



(61) 



f — 



3~& 



2G 


K n M 


5c 5 


5 


2G 


K n M 


5c 5 


5 


2G 


KnM 


5c 5 


5 


2G 


K n M 



\l~ cos 2 a + l!^ sin 2 a — ijj-g sin 2a] ai , 
[iij sin 2 a + 1^ cos 2 a + llj sin 2a] 02 , 



5c 5 5 



I33 a 3 , 



(«f - cos2 



i(lf I ) -l| ) )sin2a) 



= 0, 



2G ^ T (5) 

5c 5 11 ' 

2G Mt(5)J 
5c 5 12 > 



(62) 
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the corresponding generalized forces for M' being obtained by replacing unprimed variables by primed 
ones. The Euler-Lagrangian equations governing the dynamics of the system read now 



tii = [Eq. COJ] 



jr-g (1^ cos 2 a + I® sin 2 a — sin 2a)oi, 
2G 

a 2 = [Eq. C3] - g^O^? sin 2 a + lg> cos 2 a + sin 2a)a 2 , 



« 3 = [Eq. CE>] 
A 



2G ( 5 ) 
5^*33 fl 3. 





a 2 \ 


-1 r 


\a 2 


01/ 





ai a 2 

ft 2 01 



A-2 



-g(2igcos2a + (^-^)sin2a) 



Oi a 2 
ai a 2 

(5) t(5) 



r, 3GM' . n 
S2 H ^ — sm 2a 



CI = 



ai 

a 2 Oi 



-i r 



-2 



Ol 02 

ai a 2 

Oi 

a 2 ai 
3 GM 



A-2 
3 GM' 



Ol CL2_ 

a 2 Oi 



, sin2a-g( I gcos2a + i( I f I ) -4| ) )sm2 a ) 



[re n (a 2 (3 cos 2 a — 1 ) + a 2 (3 sin 2 a — 1 ) — a 
2 „/ 



10 r 4 

10r 5 [ K ™( a l ~ a l) sin 2a + K n( a l " a 2 2 ) sin 2a '] ~ ^ ! I 2 



+<(af(3 cos 2 a' ~ 1) + a 2 2 (3 sin 2 a' - 1) - a 3 2 )] - ^ lg> r 

r6 3GM 
-2 — 



2G 



(5) 



(63) 



Together with the corresponding equations for M' these equations describe the evolution of the binary 
system including leading order gravitational reaction. Note that in Eqs. (|63() all contributions to the 
leading order gravitational reaction are included. The quasi-static approximation for the stellar degrees 
of freedom enters into the explicit calculation of the time derivatives of the STF mass quadrupole tensor. 



3.4 Specializing to binary systems with a\ = a^, a[ = a! 2 

In the previous section we derived the Lagrangian equations of motion for arbitrary Riemann-S binaries. 
From now on we shall impose an additional constraint on the stellar degrees of freedom, requiring a\ = 
a 2 , a[ = a' 2 . On Newtonian level the equations of motion governing the dynamics of those systems have 
been derived in section 2. Now let us consider the radiation reaction part of the equations of motion. In 
principle we follow the same strategy as in the previous section. The velocity of a fluid element of star 
M relative to the star's center reads now 

u = -/3i 2 ei + flx 1 e 2 + — (xiei + x 2 e 2 ) + — x 3 e 3 , 
oi a 3 
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and the orbital motion is \i or b = —Tm^I — fM^e^- Calculating Wm and W M according to Eq. (|57Jl . the 
gravitational energy dissipation rate of the binary system yields 



_ 2G K n M 

5c 5 5 
2Gfi 



aiai(i n + i 2 2 ) + a 3 a 3 i m j - — [aiaiUn + % ) + G^lgs 



, C 5"(- I n +^^)- (64) 

The generalized dissipative forces take on a particular simple form: 

__2G«»M (5) (5) _ 2Gn n M (5) 
1 _ 5c 5 5 11 2 5 ~~ 5c 5 5 33 ' 
2G « n M l(5) 

3 5c 5 5 33 "3' 

Jvj = 0, 

_ 2G> (5 ) 
^ r ~ 5c 5 n ' 

So the Lagrangian equations of motion including leading order radiation reaction terms are given by 

.. GM' - 2 /5fciF c \ 1 3GM a x Ai G (5) 

^--ajS-^+M + -J - - Kn(1 _ |5 up + ^33 oi> 

GM' (^hLE£\l. 3GM a 3 A 3 2G (5) 
a3 ~~^ a3+ {wJJ ^~ «„(l-f) 2i?3 -5^*33 fl 3> 

01 



12 G/W 3G.M r . 2 2 \ , 1/12 /2M 2 G t(5) 
6 - — - M«l - °3) + K n("l - °3 )J " ^5 % ' 



( 5 ) 

At this point we should emphasize that in the quasi-static approximation the stellar contribution to ii-J 
vanishes, as can be easily seen from Eq. H51(l . This means, in that approximation only the time varying 
orbital mass quadrupole tensor gives rise to the emission of gravitational waves. However, coupling 
the internal dynamics of the stars with the orbital dynamics both, the orbital dynamics as well as the 
gravitational waveforms are affected by the internal dynamics. 

For our purposes it is more suitable to describe the dynamics of the binary system in Hamiltonian 
formalism. In the presence of dissipative forces the well-known Hamiltonian equations are modified 
according to 



dH . dH 

"5 j Pi — K~ 

dpi dq 



Pi = -—+F qi , (67) 



where the Hamiltonian H is given by Eq. (|29() and Eq. (|39() , respectively. The generalized momenta are 
defined in the usual way according to pi = dL/dqt. The Hamiltonian equations derived in section 2.2. 
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(cf. Eqs. (0UJI-E2IO are thus altered according to 



5 p\ GMpn n n hMP c ZGM 2 ai A x 2Gn n M T(5) 

ai+2 —r- + — - — lyai 



Pa 3 = 



2n n M a\ 
GMpn 



5r 3 



• a 3 



5r 3 nai p c 5 — n R? 5c 5 5 33 

hMP c 3GM 2 a 3 A 3 2G K n M (5) 



na 3 p c 



5-n 2R 3 5c 5 5 



*33 ° 3 ' 



Pf3 = 0, 



2G/i (5) 2 
5c 5 *12 > 



/ir 3 10 r 



1^ - al) + <(af - «?)] - l^r 



5c 5 



(68) 



while Eqs. (|38[1 remain unchanged. In the quasi-static approximation it is the time varying orbital STF 

contribi 

8GM p r 



(5) (5) 

mass quadrupole tensor alone, which contributes to ±It . Explicitly, the components of I_r read 



t(5) 

11 3 r> 

r(5) _ 2GM p. 



4^ + 3 4 + l 8 4 
r /i z p z r z 



3 r 4 



( 5) = 2GM p r 
3 r 4 



"33 



r(5) 



AGM 



r p z p z r^ 



(69) 



For numerical calculations it is useful to apply the following scaling: 



GM 



-T, p r = piCPr, p ai = PCp a% , 



GM _ 

— ^~a l: r 



GMp„ 
GM „ 



P/3 = P/3 



(70) 



Introducing the equilibrium radius i?o of a nonrotating (spherical) polytrope of mass M and polytropic 
index n the terms containing P c / ' p c can be expressed as 



kiMP c _ GM 2 /Rq 
n p c (5 — n)Ro \ R 



3/r 



(71) 
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The full set of differential equations governing the dynamics of the binary system is then given by 



5 P| Kn&i Id 
~3 "T" r ~q + 



Pax 



Pax 



Pa 3 = 



2k„Ci af 5 f 3 5-nC 2 



5 Pft K' n a[ 



1 C 2 



-I '±—L l 

2<C 2 af 5 f 3 5 - n' d 



fiifto \ 



/ _ \ 3/n 

ftp \ _ 3 ai^i 



3/n' 



_^n 03 Cl 1 
- ~3 + 



Pa, = 



5 f 3 C2 5 — n 
< ^3 , C2 1 



1 / i?n 



a[R' \R' 

3/n 



R 3 



3 a'iA[ 
R' 3 



4n n ai p r 
^Tc^f 4 



^+6^-90 



+ 



4< fflj Pr 

75 Ci f 4 



5 f 3 Ci 5 — n' 
P/3 = P^j = 



a3.R0 V ft 

1 



3 03^3 
2 i? 3 



5,oi?o \ -ft' 



2 i?' 3 



4k„ q 3 p r 
75 C 2 f 4 



4< a' 3 Pr 
75 Ci r 4 



• 6pf - 9 - 



?+6^!-9§ 



Pi 



6$ - 9 : 



p> = 


8^p0 
5 f 3 


2 

- -+ 

r 


"9^-4 




= 


1 

+ 




3 f ^~ 2 


ai = 


5 


Pax 






2K n Ci 


) 




5' 1 = 


5 


P ai 






2<C 2 


) 





16v p r 
HTf 3 



t+3p? + 18§ 



«3 = 

a 3 



Knd^ ' 
5 -., 

5 Pi 

5 

2<C 2 af ' 
r=p r , 



(72) 



— — — and C 2 = 44 = — . For the numerical calculations shown below we shall assume 



where Ci - -ttt — <um . 

1 M' n z M u 

that integration starts at periastron, i.e. the initial values for the orbit are given by 



f(0) = d(0)(l - e(0)), 0(0) = 0, PV(0) = 0, ^(0) = ^/d(0)(l-e(0) 2 ) , 

<i being the semi- major axis of the orbit. The quasi-static approximation requires all velocities inside 
the stars to be small. In particular, the mean radius ft varies only slowly with time. Hence we impose 
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i?(0) = 0, i.e. 2ai (0)5i(0)a3(0) + a!(0) 2 a3(0) = 0. This yields immediately a relation between the 
generalized momenta corresponding to ai and a 3 respectively, 



O3(0) 
'^TTTTPa! 0), 



■^Mp' (0) 
ai(0) PaiW 



Given suitable (i.e. small) values of 3j(0), 3^(0) and $(Q) respectively $'(0), we still have to fix the radii 
Rq an i?Q, respectively. Rq represents the equilibrium radius of a spherically symmetric (i.e. non- rotating) 
polytrope of mass M and polytropic index n. Assuming that the rotating star is in equilibrium at t = 0, 
we can apply the equilibrium relation between R and i?o, which has been derived in |19| 3 : 



3arcsine(0) 



e(0) 



(l-e(0) 2 ) 1/6 1 



1 



Ro = R{0) 

where e(0) = yl — (a3(0)/ai(0)) 2 is the eccentricity of the ellipsoid 



e(0) 2 e(0)arcsine(0) 



V(3- 



(73) 



4 Gravitational waveforms 



In suitable coordinates the gravitational field, observed in an asymptotically flat space, can be expressed 



as 



20 



^rad 



G 



1=2 m=- 



1-2 



D , 



^l lm (t- — )T t 



E2,lm 



(e,$) + ( - 



1-1 



«)S lm (t--)T* 2 > lm (®,$) 
c J 



(74) 

where D is the source-observer distance and the indices i, j refer to Cartesian coordinates in the asymp- 
totic space. T lm and S lm are the spherical radiative mass and current multipole moments, respectively, 
while x^ 2,lm and x^ 2 ' lm represent the so called pure-spin tensor harmonics of electric and magnetic 
type, respectively. Finally, the upper index I denotes the number of time derivatives . In leading order 
approximation only the I = 2 terms contribute to the gravitational field, which reflects the quadrupole 
character of the gravitational radiation, i.e. 



n rad 



G 



£ i 2 -(t-|)T| 2 ' 2m (e,$) 

m=-2 



(75) 



However, one might chose to use the more familiar STF mass and current multipole moments rather then 
calculating the spherical radiative ones which appear in Eq. J7SJ). In particular, the STF mass multipole 
moments are related to the spherical radiative ones by 



T m {t) 



16tt 



(2Z + 1)!! 



+ + 2) 



21(1-1) 



1/2 



1 v * lm 



(76) 



3 Of course this does not mean that the star is in equilibrium once the integration has started. 

4 Exploiting the transverse traceless character of the gravitational radiation it is possible to introduce two polarization 
vectors p and q in the plane orthogonal to the direction of propagation. This leads to the more familiar h+, h x notation, 



h 4 



PiPj 



' lilj ur ad 
n ij > 



PiQi 



where /i+ and h x represent the two polarization states. 
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where Ai = i\ . . . ii is a multi-index, while the Y l J^ are defined as 



:= (-i r( 2i - l)!! ^ (i _ 2 ^ + ro) , (4 + " • • (C + «t)C +1 

The brackets (...) denote symmetrization. Inserting this into Eq. (|76[) the spherical radiative mass 
quadrupole moments contributing to the leading order gravitational waveform read 

(77) 

-2^1^ -<£*>), (78) 
/f^-ilf-M^). (79) 

In our model only the I 20 and Z" 22 -components are present in the leading order gravitational wave field, 
since X 21 vanishes due to the symmetry of the binary system. To calculate the gravitational waveforms 
explicitly we need to know the second time derivatives of the STF mass quadrupole moment in the inertial 
frame. As before we shall compute the contributions of the star's quadrupole moments and the orbital 
terms separately. 



I 21 = 

j22 = 



Orbital contribution 

Consider inertial center of mass system. In this system the nonvanishing components of the orbital STF 
mass quadrupole tensor are given by Eq. J53J). Taking the second time derivatives and inserting the 
Newtonian equations of motion 124|) and (|25|l we end up with 



I 



(IS) 

ii 



-2/iwT0sin20 + -(1 + 3 cos 20) 
3 

GMfi 



GM 



^-(l-3cos20)0 2 



10r 3 



(1 + 3 cos 20) [n n {a\ - a 2 ) + <(af - a' 3 2 )] 



_ 2 Mi S i n 20+ —(1 + 3 cos 20) 



fir 
10r 3 



3/i 



GMfi 2 



+ ^(1-3 008 20) 



(1 + 3 cos 20) [K n {a\ - a 2 ) + <(af - af)} 



I 1 , 



r 2 



GM 



2/j,rf(j) sin 20 + — ( 1 — 3 cos 20) 

- ~ 3cos2 ^) M fl2 - a l) + <K 2 - 4 2 )] 

2 VtV^_ 20 + — — (1 — 3 cos 20) \ P 2 r - 
fir r 



^-(l + 3cos20)0 2 



3/ir : 



(1 + 3cos2</ 



GM ^l _ 3 C o S 20) [ Kn (aj - a 2 ) + <(af - a' 2 )] 



10r 3 



(80) 



(81) 
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-(IS) _ _, 

Hz — - 



2ft 
3 



GM 



2 

3// 



GMfi 2 



2 12 

r '(p 



GMfj, 
5r 3 

GM/i 
or 13 



Ii2 ^ = 2/xrr0cos20- 



2 - cos 20 + 



•2 GM 2 - 2 

r r <z> — 



[/c„(a? - 4) + <(of - a 3 2 )] 

r /■ 2 2\ i / / 12 I2\~\ 

[n n ( ai -a 3 ) +/c n (o 1 -o 3 )J , 
3GM 



(82) 



10r 3 



(K n (a 2 -a 2 )+<(af-a 3 2 )) 



/i sin 20 



p 2 . _ OM/i _ _p|_ _ 
\x r fir 2 10r 3 



(K n (af-^) + <(af-a 3 2 )) 



sin 2(f) 



(83) 



Note that tidal coupling introduces an additional contribution to I- . Since this contribution is very 
small it can be neglected for elliptical orbits and for X 22 . However, it must be taken into account when 
considering X 20 for circular orbits, where I 20 would vanish identically for a point-particle system 5 . 



Stellar contributions 

To calculate the contribution of, say, star M, let us assume the origin of the inertial system to coincide 
with the center of the star. The coordinates Xi in the inertial system and coordinates xi in the body-fixed 
system are then related by the 0(3)-transformation given in Eq. I|45l) . and the elements of the star's STF 
mass quadrupole tensor read 



(IS),star 



i s a p r refers to the body -fixed system, where the STF mass quadrupole tensor takes a particularly simple 
form. For the most general case when all semimajor axes of the star are different, the transformation to 
the inertial system yields 6 



{IS) _ 

11 — 




I22) cos 27 + 


7(^11 + ^22 — 2J 33 ), 
6 


(IS) 
22 — 




- -Z22) cos 27 


+" 7(^11 + ^22 — 2/33), 
6 


(IS) _ 
33 — 




+ ^22 — 2/33) 




(IS) _ 
12 — 




7 2 2)sin27 . 





(84) 



5 Due to tidal interaction effects the orbit is modified which in turn leads to a small, but nonvanishing X 20 -component 
of the gravitational radiation field. 

6 From now on we submit the index star. 



20 



It is easy to see that this reduces to 7-^ 
calculation yields 



(IS) _ MS) _ 1 MS) _ 1 



2^33 



(In— J33) for a\ = a 2 . A straightforward 



■MS) _ 
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cos 27 . 



In the quasi-static approximation the 2nd time derivative of the components of the stellar mass quadrupole 
tensor in the body fixed system vanishes (7^ ss 0) and we are left with 



-1 



(IS) 



MS) 

7^ S) = -2(7 n - I 22 )fl 2 sin 2 7 + [2(j u - i 22 )fl + (7 n - 7 22 )f2 

(7S) 



-2(7n - 7 22 )ft 2 cos 2 7 - 2(7n - i 22 )n + (7n - 7 22 )f2 



cos 27, 



sin 27, 



(85) 



where = 7. We should emphasize that Eqs. (|85|l imply that l\- vanishes for axially symmetric 
polytropic stars, i.e. for ai = a 2 . In other words, in the quasi-static approximation there is no direct 
contribution of the stars to the components of the gravitational held. Nevertheless the internal degrees 
of freedom give rise to modifications of the binary's gravitational wave forms due to tidally driven modi- 
fications of the orbital motion and hence of T 2m . 



Gravitational waveforms of the binary system 

As it has been pointed out before, in the case of axially symmetric stellar components of the binary 
only the orbital mass quadrupole tensor contributes to the components of hij, the actual gravitational 
waveforms being modified due to tidal coupling. Inserting Eqs. H80|) into the expressions for T 2m given 
by Eqs. (|77|) - (|79|l . the components of the leading order gravitational wave field read 

.-,2 



Z 20 = 4 



t22\ 



%(T 



3(Z 22 ) = 



'8tt 2 

-flC 



Pi 



A VrV4> ■ , 

-4 — sin 20 
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(86) 



(K„( a l ~ a 3)+ K n( a l 



a 3 2 )) > COS 20 



(87) 



(^(ai a 2 ) + <(a? fi 3 2 )) sin20 



Note the presence of an additional quadrupole-quadrupole coupling term in above equations. This term 
can be easily neglected for the 22-component of the gravitational radiation, but it will significantly 
modify the 20-component. This is due to the symmetry of the binary, the spins of the stars being aligned 
perpendicular to the orbital plane. 
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Figure 2: Example for the stellar oscillation within the affine model. Plotted is the oscillation of the 
semi-major axis a± of star 1 in a tidally coupled binary system. 



5 Numerical results and discussion 

In the last few years there has been a growing interest in investigating the influence of tidal interaction 
onto the inspiral process and the actual leading order gravitational waveforms emitted by a binary 
system. Within the framework of linear perturbation theory of stellar oscillations, close binary systems 
of nonrotating neutron stars where studied by Kokkotas and Schafer JJj and later on by Lai and Ho [Hj • In 
these papers it was shown that tidal interaction may draw energy from the orbital motion, thus speeding 
up the inspiral process. This effect is strongest in the case of a tidal resonance when the orbital frequency 
is the mth fraction of a star's eigenfrequency, m being an integer number. Resonant tidal oscillations of a 
nonrotating white dwarf-compact object binary were investigated by Rathore et al. |21| . but in this paper 
dissipation due to gravitational waves emission was not included. On the contrary our model incorporates 
not only dissipation due to gravitational wave emission, but also stellar rotation. Moreover, we do not 
restrict ourselves to small, linear perturbations of thermal equilibrium, but allow for oscillations with 
arbitrary large amplitude. However, there is only one oscillation mode per star present in our analysis 
(this mode corresponding roughly to the f-mode), this means we truncate the all the higher eigenmodes 
which are incorporated in the linear perturbation scenarios. If much of the oscillation energy is stored in 
these higher eigenmodes the gravitational waveforms calculated from our model will fail to give results 
of high accuracy. Therefore we restrict ourselves to moderate orbital separations, when the dominant 
contribution of the oscillation energy comes from the /-mode. In Fig. |2]the oscillation of the semi-major 
axis ai is shown for a particular example. 

The most prominent feature of the model is that it allows us to study the long term evolution of the binary 
system. In fact, in our numerical calculations we were able to follow the orbital evolution over hundreds 
of periods. In particular we studied the influence of the equation of state on the binary's dynamics and 
thus onto the leading order gravitational wave pattern emitted during the time evolution. In the past it 
was argued that the influence of the polytropic index will not be reflected in the gravitational waveforms 
almost until the final plunge down [3]. While this is certainly true if one considers only a few orbits 
this suggestion has to be carefully checked for long term evolutions. Given suitable initial conditions 
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our model allows for the investigation of the orbital evolution over hundreds of periods. It is thus an 
ideal tool to study the influence of the polytropic index on the dynamics as well as on the gravitational 
waveform in great detail. Note that the equations of motion given in Eqs. (|72[1 do not apply to the limit 
of incompressible fluids (n — 0) and to the relativistic limit, which is represented by T — 4/3, i.e. n = 3. 
For n = Eq. I|71|l has to be replaced by 

5fci P c P c 
— 1_£ = 2—. 

nn n Pc Pc 

(For an explicit expression of P c /p c in terms of the stellar degrees of freedom see 22).) For n — > 3 
numerical integration becomes increasingly instable. Moreover, Eq. I|73|l . which gives a relation between 
the equilibrium radius i?o of a nonrotating polytrope and the equilibrium mean radius R of a rotating 
polytrope, becomes singular for n — 3. 

In Fig. [3]we compare the orbital evolution of a slightly elliptic (e(0) = 0.4) equal mass binary system for 
different choices of n' . As expected tidal coupling induces a periastron advance, but it is clearly visible 
that the periastron advance increases with decreasing values of n'. Not too surprisingly, the polytropic 
index n' also affects the inspiral process. For fixed value of n the inspiral process speeds up the smaller 
the values of n', i.e. the larger the value of V . As it is shown in Figs. ^JandElthe effect of the equation 
of state on the gravitational waveforms emitted by the binary cannot be neglected even for moderate 
orbital distances if one considers the long-term evolution of the binary system. For circular orbits, the 
tidally induced modification of the gravitational wave pattern is strongest in the I 20 — component of hij 
(see Fig. |5jl, but also in the 22— component there is an significant phase shift due to different polytropic 
indizes n! . This is demonstrated in Fig. 0] where we compare the 5R(X 22 )— component of the gravitaional 
wave field for different values of n'. In these particular examples the orbit is assumed to be circular in 
the absence of tidal perturbation and the initial orbital separation is taken to be d(0) — 50. Even in this 
case a phase shift in the 5R(X 22 )— component is already obvious after 10 orbital periods. 
As already mentioned before, for a circular orbit the influence of the equation of state is strongest reflected 
by the the X 20 -component of hij. For a point particle binary, where the leading order gravitational 
waveform is known analytically, this contribution to hij would vanish for a circular orbit. This is not the 
case if the orbital motion is tidally coupled to the internal motion of the stars (see Fig. 03 . 
For elliptic binaries the influence of the polytropic index n' on the I 22 -component is demonstrated in Figs. 
Elandd Note that even for large orbital distances different polytropic indices n' manifest themselves in a 
remarkable phase shift of the X 22 -component of h^ already after a few orbital cycles. This is demonstrated 
in Fig. [7| where we compare 3?(X 22 ) for different values of n' at different stages of the inspiral process. 
To summarize, we have shown that including the internal structure of the stars will lead to significant 
changes of the leading order gravitational wave pattern compared to a point particle binary. Moreover, 
even for moderate relative distances of the stars there is a notable phase shift in the X 22 -component of the 
gravitational radiation field after some orbital periods. Of course our model does not respect all features of 
stellar dynamics. In particular it does not account for higher order oscillation modes which are present in 
linear perturbation theory. Nevertheless, our model could give an approximative description of a neutron 
star-neutron star binary. Although the internal structure of a neutron star is still not understood very 
good (the best known equations of motion are only given in a tabulated form), astrophysical observations 
indicate that the thermodynamics of neutron stars can be approximately described by a polytropic 
equation of state with polytropic index n ~ 0.5 — 1.0 (see e.g. (2U)- Moreover, all low-mass white dwarfs 
can be modelled with an effective polytropic index n « 1.5. Thus our model might be applicable to study 
the long term evolution of NS-NS binaries, NS- White Dwarf binaries or binaries consisting of either a 
neutron star or a White Dwarf and a compact object. 
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A The Lane-Emden equation reviewed 

For spherically symmetric polytropes des equations of hydrodynamics are given by 

Gm{r)p{r) 



dP 

dr r 2 
which can be combined to a single equation 

1 d fr 2 dP 
r 2 dr \ p dr 

Introducing dimcnsionlcss variables £ and 6 as 

p:=PcB n , 

with 



^ m A I \ 2 

— — = VKp\r)r 
dr 



= -iirGp. 



a 



[n + l)KpS~ 



1/2 



AttG 



we obtain the well known Lane-Emden equation 

,d8^ 



1 d 



0(0) - 1, 



9'(0) = 0. 



(89) 



(90) 



For n < 5 the solution of Eq. H9()|) decreases monotonically and becomes zero at a finite value of £. This 
value, denoted by £i, is characterized by vanishing pressure and density and hence represents the star's 
surface. 



B Chandrasekhar's coefficients Ai and J 

The quantity J is defined in JHj (chapter 3) according to 



f°° du 

J := a ia2 a 3 / — = (91) 

Jo y / {al+u)(a2 + u){a 3 + u) 



while the coefficients Ai are given by 



A- [°° ^ U (92) 

Jo (a 2 + u) ^ (aj + u) (a 2 , + u) (a| + u) 



One can easily show that 

J = a\Ai + a\A 2 + alA 3 . 
Another useful relation that can be easily derived from Eq. I|91ll is 

dJ 1 ,„ 2 



\J-a\Ai) (no sum) (93) 
0^ ai 



25 



and 

87 1 87 1 

= —{J - a\M), -i~ = ~{J ~ 4A 3 ) for ai =a 2 . (94) 

ocli ai 8a 3 a 3 

In the special case a\ = a 2 > 03 it is possible to find an analytical expression for A\ and A 3 . Solving 
yields 



y/l - e 2 1 - e 2 
A\ = A2 = 5 arcsin e - — (95) 

e J er 

2 2v / r _ e 2 

A 3 = — - arcsin e (96) 



where we defined 



2 

03 ' 



e:=^/l-(^). (97) 



j(5),star |3 e y 0n( j ^ e quasi-static approximation 



b 

Throughout the paper we applied the quasi-static approximation to calculate the radiation reaction terms 
of the stellar degrees of freedom. This approximation is justified if all internal velocities and accelerations 
inside the star are small. Dropping this assumption, the result is much more complicated. In particular, 
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the 5th time derivative of the STF mass quadrupole tensor of M in the corotating coordinate frame reads 



^ = (^1-122) 



(i6^ 5 - 4on 2 n - 6onn 2 + n {4) ) 




- (80O 3 o - 2onn - iof^ (3) ) 



cos 2a 
— sin 2a 


- cos 2a sin 2a 
sin 2a cos 2a 




+ (±11 - ±22) 



cos 2a sin 2a 



(40ft 4 - 30ft 2 - 40ftft) sin 2a -cos 2a 



+ 10(1^-1 



+ 50(^-4^) 





(— sin 2a 
cos 2a 


sin 2a 
cos 2a — 


— cos 2a sin 2a 
sin 2a cos 2a 










cos 2a 
sin 2a 






27 




Figure 3: Orbital evolution and periastron advance due to tidal coupling for an elliptic, equal mass 
binary with initial values e(0) = 0.4, d(0) = 50,ai(0) = a[(0) = 10,a 3 (0) = a 3 (0) = 8,/3(0) = /?'(0) = 
0.01, ai(0) = 0.02, a[(0) = 0.01. All systems have polytropic index n = 1, while star 2 obeys a polytropic 
equation of state with n' = 0.2 (upper left), n' = 0.5 (upper right), n' = 1 (lower left) and n' = 1.5 (lower 
right). Star 1 has polytropic index n = 1 in all scenarios. 
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Figure 4: The 3?(i 22 )-component of the leading order gravitational waveform for a circular orbit (e(0) = 0) 
for polytropic indices n! = 0.2 (dotted line) and n' — 2 (solid line) at different time stages. The polytropic 
index of star 1 is assumed to be n = 1, the initial parameters are the same as in Fig. [3] 
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Figure 5: I 20 -component of the leading order gravitational waves for a binary with polytropic index 
n = 1 for star 1 and n' = 0.2 (upper left), n' = 1 (upper right), n' — 1.5 (bottom left), n' = 2 (bottom 
right). The eccentricity at t = is assumed to be zero, the initial parameters are the same as in Fig. 0] 
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Figure 6: Influence of the polytropic index nf on the real part of I 22 for an elliptic binary with e(0) = 0.4. 
The initial parameters are the same as in Fig. 0] Left: n' — 0.2, Right: nf = 1.5. 




Figure 7: 5R(X 22 )-component of the leading order gravitational radiation field emitted by an elliptic, 
equal-mass binary (e(0) = 0.4). Shown are the waveforms for n' = 1 (solid line) and n' = 0.2 (dotted 
line), the polytropic index of star 1 is assmued to be n = 1. The semi- major axis at the beginning of 
the orbital evolution is taken to be d(0) ~ 100, while the initial parameters for the stars are given by 
oi(0) = ai(0) = 20, a 3 (0) = a' 3 (0) = 19, ai(0) = a'^0) = and /3(0) = 0.002, /3'(0) = -0.002. 
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